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Abstract 

We present a method of analyzing invariant-mass spectra of kaon pairs re¬ 
sulting from decay of 4> mesons produced in high-energy heavy-ion collisions. 
It can be used to extract the shifts in the mass and the width (AM and AF) 
of the (j) mesons when they are inside the dense matter formed in these colli¬ 
sions. We illustrate our method with the help of available preliminary data. 
Extracted values of AM and AF are significantly larger than those obtained 
with an earlier method. Our results are consistent with the experimentally 
observed px dependence of the mass shift. Finally, we present a phenomeno¬ 
logical relation between AM and AF. It provides a useful constraint on 
theories which predict the values of these two quantities. 
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1 Introduction 


Properties such as mass, decay width, size, etc. of hadrons in general and light vec¬ 
tor mesons (p, cn, 0) in particular, are thought to undergo a change when the surround¬ 
ing medium is raised to a higher temperature and/or density; for recent reviews of hnite- 
temperature and hnite-density effects, see Refs. [Q and 0, respectively. Interest in this 
scenario arises largely due to its close connection with phase transitions in quantum chromo¬ 
dynamics (QCD). Recent years have witnessed intense experimental and theoretical activity 
in this area In this paper, we focus on the mass as well as the width of the 0 meson. 

Recently E-802 experimental collaboration at AGS (BNL) has reported preliminary re¬ 
sults on the shift in the mass of the 0 mesons produced in central Si -|- Au collisions at 
14.6 A GeV/c |^J^. The mass was determined from the invariant-mass distribution of kaon 
pairs arising from the dominant decay mode 0 —> K^K~. It was studied as a function of 
the multiplicity of the collision events and the transverse momentum px of the accepted 0’s. 
It was found that in the events with the highest multiplicity (top 2% target-multiplicity- 
array TMA cut) the mass drops by2.3±0.9±0.1 MeV compared to the free-space value 
1019.413 ± 0.008 MeV [^. For the 0 mesons with pt < 1.25 GeV/c, the shift was even 
more (3.3 ± 1.0 ± 0.1 MeV), whereas for p^ > 1.25 GeV/c, there was no apparent shift. No 
numbers were reported for the shift in the decay width of the 0 mesons. However, from the 
conhdence contours for 1, 2 and 3 standard deviations in the observed mass versus width 
plot given in P] , it appears that the central value of the width is higher than the free-space 
value 4.43 ± 0.06 MeV 0 by about 0.78 MeV. If confirmed, this will be the first evidence of 
the modification of hadronic properties inside dense matter formed in heavy-ion collisions. 

In a more recent publication [0, the E-802 collaboration has reported no shift in the mass 
and the width of the 0 meson. This, however, does not contradict the earlier observations 
P because in Ref. []^ the top 7% of the highest-multiplicity events were clubbed together 
and no multiplicity dependence of the mass was reported. In the earlier work [Q, the shift 
in the mass was seen when the TMA cut was applied at 2%. 
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2 Method 


In this paper, we present a different method of analyzing the invariant-mass spectra of 
kaon pairs. We then use the data in Refs. [4-5] only to illustrate our method. We find 
that our method yields shifts in the 0-meson mass and width which are several times larger 
than the values quoted above. We caution the reader that since these data are preliminary, 
our numerical results are obviously subject to change. However, our main point is the 
method presented here which is independent of whether the published data [4-5] are eventually 
confirmed or not. We think the present method has relevance to the analysis of future data 
on 0 production. This acquires added importance in view of the possibility of a similar 
experiment being performed at CERN APS [§. 

Let Mq and Mi be the rest masses and let Pq and Pi be the widths of the 0 meson in 
the free space and in the dense medium, respectively. We define the shifts as 

AM = Mo-Ml, AP = Pi-Po, (1) 

so that both are positive when the mass drops and the width increases with respect to their 
free-space values. 

In Ref. [^, the invariant-mass spectrum of the K^K~ pairs was fitted by a function con¬ 
sisting of a background term and a relativistic Breit-Wigner (BW) resonant term convoluted 
with a Gaussian experimental mass resolution function. This procedure yielded the values 
of the shifts AM and AT given above. A background-subtracted mass spectrum is displayed 
in Fig. 1 where the histogram corresponds to 3% TMA and low-p^ cuts and is taken from 
Fig. 2(b) of Ref. 0. The dashed and solid curves represent single BW resonance terms 
convoluted with a Gaussian as stated above. Areas under the two curves are normalized 
to that under the histogram. The dashed curve corresponds to the scenario where there is 
no shift in the mass and the width. The solid curve corresponds to AM = 2.3 MeV and 
AT = 0.78 MeV. Glearly the data indicate a shift in the mass and the width of the 0 meson. 
However, are the values of AM and AT obtained by fitting a single BW resonance term to 
the (background-subtracted) data correct? We think they are not. 
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Since the mean lifetime of 0 in its rest frame is about 45 fm/c, a majority of 0’s are 
expected to decay long after the dense medium in which they were produced has ceased to 
exist. That is, they will decay essentially in free space (AM = 0 = AT). The rest of the 
0’s, however, decay while still inside the dense medium. Hence a better procedure would be 
to fit the background-subtracted data with two instead of one BW terms, one unshifted and 
another shifted, added with appropriate weights. This we now proceed to do. 

We work in the center-of-mass frame of the dense medium formed in the nucleus-nucleus 
collision. We employ the natural units c = h = 1. Let / be the fraction of 0’s decaying 
inside the medium; then (1 — /) is the fraction decaying in free space. We reanalyze the 
background-subtracted data of Ref. by htting them with two BW terms — an unshifted 
BW with mass Mq, width Tq and weight (1 — /) and a shifted BW with mass Mi, width Ti 
and weight /: 


dNx+K- 

dM 


(1 - /) BW,iM, Mo, To) + / BW,{M, Mi, Ti). 


( 2 ) 


Here BWc denotes the relativistic Breit-Wigner resonance term convoluted with a Gaussian 
experimental mass resolution function: 


RW,(M, Mo, To) = J BW{M', Mo, To) exp 

Mo ro(M) 2M 


1 (M -M'^ 


a 


dM' 

■ 71 ^’ 


a 


RW(M,Mo,ro) = 


( 3 ) 

( 4 ) 


TT (M2 - rg(M) ■ 

The experimental mass resolution a is taken to be 2.2 MeV |0. The energy-dependent width 
ro(M) is taken to be [|T^ 


ro(M) = To ^ 




r + ?o 

where the kaon momenta q and go are given by 


( 5 ) 


q = {M''-li- Mlyl^ and 9o = (V?/4 - 


( 6 ) 


Expressions for BWc{M, Mi, Ti), BW{M, Mi, Ti) and ri(M) are obtained by replacing Mo 
by Ml and To by Ti in Eqs. (3-6). We hnd that the numerical results with energy-dependent 
and energy-independent widths to be practically the same. 
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We now derive an expression for /. Let tq = Lg ^ and ri = be the mean lifetimes of 
0 at rest in the free space and in the dense medium, respectively. Then the corresponding 
lifetimes when the 0 is in motion would be ro 7 and ri 7 , due to the time dilation. Here 
7 = (1 — is the usual Lorentz factor and v is the average velocity of the 0’s in the 

medium. Let d be the “size” or “radius” or a typical linear dimension associated with the 
extent of the medium. The average time the 0’s produced in the medium would take to 
traverse this distance is d/v. If Nq is the number of 0’s at time f = 0, then at time t = d/v, 
only iVo exp(—Tid/ny) of them would be left. Hence / = 1 — exp(—Tid/ny). Now 

ny = p/Mi = \Je‘^ - Mf/Mi = \J M|. cosh^i/cm - /Mi, 

where p, E, Mt and Ucm denote, respectively, the momentum, energy, transverse mass and 
center-of-mass rapidity of the 0 traversing the medium. Hence 

/(AM, AT) = 1 - exp(-Mirid/7M/ cosh^p,,, - Mf). (7) 

Interactions of the outgoing with the medium were ignored in Refs. [4-5] and we 
too shall ignore them. This is a good hrst approximation because if these interactions were 
important the 0 peak would have been washed out, whereas experimentally, a distinct, 
narrow peak is seen. Given the quality of the available data, we think, the present approach 
is adequate. In view of the possibility of better-quality data becoming available in the future, 
it is desirable to include these interactions by performing a detailed Monte-Carlo calculation. 
This would entail a considerable amount of extra work, and we leave it to the future. Such 
an elaborate calculation may change the values of AM and AT to some extent, but we do 
not expect it to change our conclusion that these will be more realistic than those resulting 
from the one-BW £t to the data. If the kaon-medium interaction is ignored, then Mt and 
i/cm of the 0 are the same as those of the pair. In a Monte-Carlo simulation one can 
extract them from experimental data by modeling the interaction. 
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3 Results and discussion 


Results of a least-squares fit to the same experimental data as in Fig. 1, with two instead 
of one BW terms, are shown in Fig. 2. The two dashed curves in Fig. 2 correspond to the 
two convoluted BW terms, one shifted and the other unshifted. The solid curve corresponds 
to their weighted sum as in Eqs. (2-6). The input parameters are d = 5 fm, Mt = 1100 MeV 
and Hem = 0.3 [Q. The two htted parameters are AM = 6.0 ± 1.9 MeV and AT = 5.4 ± 3.9 
MeV. The resultant value of / is 0.37. If Mt = 1200 MeV, the two htted parameters are 
AM = 6.2 ± 2.4 MeV and AT = 7.7 ± 5.0 MeV. The resultant value of / is 0.34. Note that 
the above values of the shifts are several times larger than those obtained with a one-BW 
ht. Secondly, even for an unrealistically large value of d, say 10 fm, the fraction / of 0’s 
decaying inside the dense matter is ~ 0.4 to 0.5 showing that the assumption of a single 
BW resonance is questionable. 

It is evident from Eq. (7) that as pt increases the fraction / decreases. This means a 
larger fraction of 0’s decay in free space and hence the mass shift decreases. This is exactly 
what has been observed experimentally [^. 

Finally, we present an empirical relation between AM and AF. We dehne 

M = (1 - /)Mo + /Ml, (8) 

AM = Mo - M. (9) 

On substituting Eqs. (1) and (9) in Eq. (8) and simplifying we get 

AM = /(AM, AF)AM. (10) 

Given the experimental histogram as in Fig. 1, it is straightforward to determine its centroid 
and it is reasonable to approximate M by the centroid. Equation (9) can then be used to 
determine the “apparent” shift AM. For a hxed AM, Eq. (10) yields a curve AF versus 
AM. We illustrate this in Fig. 3 for AM = 1, 2, 3 and 4 MeV. The input parameters are 
as above, namely d = 5 fm, Mt = 1200 MeV and Hcm = 0.3. Interestingly, as AM increases, 
AF decreases. This is easy to understand from Eq. (10) because for a hxed AM, as AM 
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increases, / has to decrease, which from Eq. (7) requires Ei and hence AE to decrease. 
The physics of this is also clear if one considers the three curves in Fig. 2. Equation (10) 
together with the experimentally determined AM provides a useful constraint on theories 
which predict the values of AM and AE. For example, if the apparent shift AM extracted 
from experimental data is 2 MeV, then any model which claims to explain the data should 
have AM and AE not inconsistent with the curve labeled “2” in Fig. 3. 

A variety of approaches, namely lattice QCD, QCD sum rules, Nambu-Jona-Lasinio 
model, quantum hadrodynamics, bag models, instanton-liquid model, etc. have been used 
in the literature to predict the masses and widths of the light vector mesons at hnite tem¬ 
peratures (T) and/or densities (p). Here we focus on some of the recent calculations for the 
0 meson. Asakawa and Ko [|^ have used the QCD sum rules to calculate the 0-meson mass 
when both T and p are hnite. For the values of T and p considered by them, the predicted 
drop in the mass is as high as a few hundred MeV. They have proposed a double-peak 
structure in the invariant-mass spectra of dilepton pairs arising from 0 decay as a signal of 
the transition from quark-gluon plasma to hadronic matter. Finite-temperature effects on 
the 0-meson mass and width have been studied recently by Shuryak and Thorsson, Haglin 
and Gale, Bhattacharyya et ah and Song . 


4 Conclusions 

In conclusion, (z) we have presented a method of analyzing invariant-mass spectrum of the 
kaon pairs; the resultant values of the shifts in the 0-meson mass and width are signihcantly 
larger than those obtained with an earlier method. This difference arises due to the long 
lifetime of the 0 which ensures that only a small fraction of them decay inside the dense 
matter, [ii) The model presented here is consistent with the experimentally observed pr 
dependence of the mass shift, {iii) We have presented a phenomenological relation between 
the shifts in the mass and the width. 

It is worth reiterating that the data in [Q are preliminary. In view of the interesting 
conclusions that could be drawn from the data, it would be useful to have a more thorough 
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experimental investigation of the shift in the mass as well as the width of the 0 meson. One 
would like to have data with better statistics, on a variety of targets and beams, and at 
various energies. One would like to see invariant-mass spectra of not only the kaon pairs 
but also the lepton pairs resulting from 0 decay. 
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FIGURES 

FIG. 1. The background-subtracted invariant-mass spectrum for the K~^K~ pairs. The his¬ 
togram represents preliminary experimental data from [5]. The dashed and solid curves represent 
single BW resonance terms convoluted with a Gaussian as explained in the text. The dashed curve 
corresponds to the situation where there is no shift in the i;i)-meson mass and width. The solid 
curve corresponds to AM = 2.3 MeV and AT = 0.78 MeV. 

FIG. 2. The background-subtracted invariant-mass spectrum for the K^K~ pairs. Experimen¬ 
tal data as in Fig. 1. The two dashed curves correspond to the two BW terms, one unshifted and 
the other shifted. The solid curve corresponds to their sum weighted by the factors (1 — /) and /, 
respectively; see Eqs. (2-6). 

FIG. 3. For a fixed AM, Eq. (10) provides a relation between AF and AM. This relation is 
represented by the solid lines labeled by the value of AM in MeV. 
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